Hydrodynamic bubble coarsening in off-critical vapour-liquid phase separation 
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Late-stage coarsening in off-critical vapour-liquid phase separation is re-examined. In the limit of 
bubbles of vapour distributed throughout a continuous liquid phase, it is argued that coarsening 
proceeds via inertial hydrodynamic bubble collapse. This replaces the Lifshitz-Slyozov-Wagner 
mechanism seen in binary liquid mixtures. The arguments are strongly supported by simulations in 
two dimensions using a novel single-component soft sphere fluid. 
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Late stage coarsening kinetics in phase separating 
mixtures shows a wealth of interesting non-linear cou- 
pled phenomena |p]-|lo|]. If dynamical scaling holds, the 
growth may be characterised by a scaling law relating 
the typical size of domains, L, to the elapsed time t, 
thus L ^ t". The standard picture of coarsening in a 
space dimension d = 3 critical quench (equal propor- 
tions of the two phases) is as follows [|l| . An initial diffu- 
sive regime (a — 1/3) is followed by accelerated growth 
when hydrodynamic modes become activated. If a fluid 
is characterised by its density p, viscosity rj, and surface 
tension a, then hydrodynamically limited coarsening fol- 
lows a law like L/Lh = f{t/TH) on dimensional grounds, 
where Lh = r]^ / per and Th = rf / pa^ are hydrodynamic 
length and time scales. If scaling holds, f{x) ~ x". 
The expectation is that an initial viscous hydrodynamic 
(VH) regime {a — \) crosses over to an inertial hydrody- 
namic (IH) regime (a — 2/3), since the Reynolds num- 



ber Re = pLL/ri = //' 



^ grows like x in the VH 



regime and eventually inertial effects must become im- 
portant. The observation that Re continues to grow in 
the IH regime (as x^/^) has recently prompted Grant 
and Elder to argue that turbulent remixing will limit the 
asymptotic exponent to a < 1/2 However the cur- 
rently available numerical evidence suggests a breakdown 
of pure scaling is more likely, in the sense that coarsening 
is no longer dominated by a single length scale . 

In an off-critical quench (unequal proportions of the 
two phases) the situation seems much clearer espe- 
cially if the volume fraction of the minority phase is van- 
ishingly small j^. For the case of a binary liquid mix- 
ture, one rapidly establishes a population of droplets of 
the minority phase distributed throughout the continu- 
ous phase. Such droplets can coarsen via the Lifshitz- 
Slyozov-Wagner (LSW) mechanism in which the larger 
drops grow at the expense of the smaller drops, as mate- 
rial diffuses through the continuous phase driven by weak 
chemical potential gradients. Diffusion limited droplet 
coalescence may additionally occur. Both mechanisms 
lead to a = 1/3 in d = 3. 

Up to now I have not drawn a clear distinction between 
coarsening in a single component system quenched into 



vapour-liquid coexistence, and demixing in a binary liq- 
uid mixture. Both have similar conserved order param- 
eters, and one expects similar kinetic coarsening mecha- 
nisms to hold for critical quenches. In this Letter though, 
I argue that coarsening in vapour-liquid systems is funda- 
mentally different from binary liquid mixtures in the limit 
of bubbles of vapour distributed throughout a continu- 
ous liquid phase. This reflects a fundamental difference 
in the way the conserved order parameter is transported, 
and is motivated by simulations of coarsening in a novel 
single-component soft sphere fluid described below. 

Visual observation of the bubble coarsening regime 
shows that an LSW-like mechanism operates, but much 
faster than would be expected from diffusion of the order 
parameter. The large bubbles grow by bulk movement 
of fluid, at the expense of the small bubbles which col- 
lapse under the influence of their surface tension. This 
is strongly suggestive of a new hydrodynamically limited 
bubble coarsening regime, in which flow fields are gen- 
erated by pressure gradients caused by a distribution of 
bubble sizes. These pressure gradients are analogous to 
the chemical potential gradients in the LSW mechanism, 
but in this case the continuous fluid phase responds hy- 
drodynamically. 

To gain insight, I consider the paradigmatic problem of 
an isolated bubble, initial radius Rq, which collapses from 
rest under the influence of its surface tension, in zero am- 
bient pressure. Related problems of bubble growth and 
collapse under varying ambient pressure fields go back 
to Rayleigh Q , and are of great technological and mil- 
itary significance in cavitation phenomena |l^ ]. In these 
problems, the flow field is radial, v = {R^R/r^)r, where 
R{t) is the drop radius. Incompressibility requires the 
1/r^ dependence, and the fluid velocity is matched to 
the bubble radius velocity a,t r = R. This flow field is 
irrotational and can be derived from a velocity potential 
(f) = —R^R/r. The drop radius obeys the equation 



p{R + 3if/2) = p{R) - p{oo) 



(1) 



where p(r) is the pressure. This can be derived for ex- 
ample by applying Bernoulli's principle to a streamline 
extending from the surface of the drop to infinity p3|. 
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In our case p{R) = —2a/R from the Laplace law (the 
pressure inside the droplet is « 0), and p{oo) = 0. The 
above equation can be integrated to find that the bubble 
collapses in a time 



^collaps 



dx 



L)l/2 



(2) 



(the integral can be evaluated in terms of Euler's B- 
function to be B(|,i)/2 = 0.874). This result can be 
rewritten as (tcoUapse/?//)^^'' c< Rq/Lh- It is not sur- 
prising this coincides with the IH scaling law, since we 
are only considering a special case of the Navier-Stokes 
equations in which viscosity effects are absent. 

This result for isolated bubble collapse suggests that in 
the LSW-like hydrodynamically limited bubble coarsen- 
ing regime, the mean bubble size grows as i/f (t/T'^r)^/'^, 
identical to the growth law for IH coarsening. What is 
perhaps surprising is the indication that there is no VH 
regime, since isolated bubble collapse is purely inertial 
(but see discussion below). The above argument is for 
d = 3, but repeating the analysis in c? = 2 obtains the 
same scaling law with additional logarithmic corrections. 

I now describe the simulations which motivated the 
above arguments. These were undertaken with a novel 
soft sphere fluid based on the dissipative particle dynam- 
ics (DPD) method Q. Two features make DPD attrac- 
tive for phase ordering studies: firstly one has a fluid 
of soft spheres which allows comparatively large time 
steps to be taken in the integration algorithm, secondly 
a momentum-conserving thermostat is used which pre- 
serves hydrodynamic modes. 

The soft spheres in DPD interact with short-range 
forces of the form fy = A{1 — rij/rc)vij, acting between 
all pairs of particles i and j for which < Tc, where 
the positions are and r^, Vij — Tj — r^, rij — |ry| 
and Tij — r.ij/rij. Such a force law leads to a predomi- 
nantly quadratic equation of state (EOS) which cannot 
be engineered conveniently to have a van der Waals loop. 
Phase ordering studies, using DPD in its original form, 
are therefore limited to binary fluid mixtures where it 
has been applied with some success 

Usually a van der Waals loop results from a hard core 
repulsion combined with a long range attraction, for ex- 
ample in previous studies using Lennard- Jones potentials 
|^,|o|. However, introducing hard cores into DPD spoils 
the attractiveness of the method. Therefore an alterna- 
tive method was pursued which may perhaps be termed 
many-body DPD. 

In many-body DPD, the amplitude of the interaction 
law is made to depend on a weighted local density, sam- 
pled at the particle positions in an additional sweep 
through the pairs of interacting particles. To be spe- 
cific, the following quantity is computed for each particle: 
Pi = J2j=i '^{'''ij) (* = 1 ■ ■ • N) where w{r) is a suitably 
normalised weight function vanishing for r > Tc- Then 



the amplitude of the interaction between a pair of parti- 
cles i and j is made to depend symmetrically on pj and . 
In this way, almost completely arbitrary density depen- 
dence can be introduced into the EOS, including a van 
der Waals loop. Many-body DPD in this form has been 
described by Pagonabarraga and Frenkel ||l5|, and was 
also invented independently by Groot p6| . Such density- 
dependent potentials have featured recently in coarse- 
grained models of polymers , and in principle can be 
constructed to have arbitrary phase behaviour ||l^ . 

The specific model used in the simulations comprises 
DPD with an additional interaction force of the form 
fcxtra ^ ^- _ {r,, < r^). The d = 2 

weight function is w;(r) = (6/7rr^)(l— r/r^)^ (r < r^) [ p^ . 
The EOS for this model is then predominantly cubic. By 
making A < 0, B > and < Tc, one can produce a 
van der Waals loop and induce vapour-liquid coexistence 
with a sharp interface of width ~ Vc (see Fig. |l|). Three 
suitable parameter sets were identified for d — 2 simu- 
lations, and are given in Table | along with the density 
and viscosity of the liquid phase, and interfacial tension 
computed by separate simulations. In common with pre- 
vious studies, I fix the units by setting m = = kT = 1 
(to is the mass of the particles). 

Armed with this model, d = 2 simulations of vapour- 
liquid phase ordering were undertaken by preparing a 
random distribution of particles in a square domain at 
a density (j)p where is the desired liquid phase area 
fraction. Typical simulations comprised 5-10 runs of 
lO'^-lO'^ particles in a domain of side lOOr-c, at = 0.9. 
These were run out to times of order 150 (DPD units) 
where only one or two large vapour bubbles remain. The 
progress of the simulation was monitored by overlaying 
a grid at a resolution 0.5rc, computing the local density 
in each grid element (pixel), and determining connected 
sets of pixels where the density is less than a critical value 
(usually p = 1) by cluster analysis (see Fig. |^). These 
clusters are identified with bubbles, and statistics on the 
bubble size distribution were obtained. Most attention 
was paid to how the mean bubble area (A) grows with 
time; very similar results were obtained with {A'^^'^). 

Fig. m shows data from runs for all three parameter 
sets, reduced using the hydrodynamic length and time 
scales. There are two key points to note. Firstly the 
data collapses onto a single curve despite the order of 
magnitude variation in Lh and Th across the sets. This 
is a strong indication that bubble coarsening is hydro- 
dynamically limited. Secondly the data follows a scaling 
law virtually indistinguishable from that expected in an 
IH regime over nearly four decades in reduced time. This 
is a strong indication that the above arguments concern- 
ing inertia-dominated bubble collapse are correct. To be 
precise, the scaling law found by fitting across all data 
sets is (A) = (4.2 ± 0.5)L|^(i/'7H)i-30±o o2. The barely 
significant deviation of the exponent from 2a = 4/3 may 
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be due to the logarithmic corrections expected in d = 2 
or the remnant effects of viscosity. 

By way of comparison I have also done some corre- 
sponding simulations on a binary liquid mixture. This 
is achieved in regular DPD by introducing two species of 
particle with an additional repulsion between them 0. 
The parameter sets in Tabl e p were used, and the re- 
sults are also shown in Fig. M. The growth exponent is 
that expected for diffusion-limited droplet coalescence in 
d = 2 as has been seen in previous work 1^. Note the 
absence of a scaling collapse between the two parameter 
sets in this case. 

It may be worthwhile to interject a brief comment 
about the opposite limit to the bubble coarsening regime, 
namely a dispersion of drops in a vapour phase. If 
the vapour is sufficiently dilute, coarsening can pro- 
cees by ballistic drop coalescence, as has been reported 
previously for a Lennard- Jones fluid ||l0[| . Following 
Bray p]|, si mple scaling arguments suggest a growth law 
[(j)t^/kT/pY^^'^^^K I have conffi'med the vahdity of var- 
ious aspects of this growth law in d = 2 at (/) = 0.1 and 
0.05. 

In summary therefore, I have presented arguments 
which show a fundamental difference in phase coarsening 
between binary liquid demixing and vapour-liquid phase 
separation. Simulations of a soft sphere fluid strongly 
support the idea of an incrtial hydrodynamic (IH) bub- 
ble coarsening regime for an off-critical vapour-liquid 
quench. The relevance of the critique of Grant and Elder 
is not immediately apparent in this case, nevertheless 
there has to be a mechanism to dissipate the interface en- 
ergy. It may be that it is converted to kinetic energy in 
the fluid, and subsequently lost in viscous dissipation by a 
turbulent cascade, as suggested for the critical quench in 
a binary liquid mixture g. But, in real cavitation, a sig- 
nificant amount of the interface energy is converted into 
pressure waves . Thus another possibility is that fluid 
compressibility plays a role, and the interface energy is 
(partly) dissipated through sound waves. The subtleties 
of these mechanisms do not show up in the statistics of 
the bubble size distribution since the data in Fig. ^ in- 
dicates the mean bubble size follows the IH growth law 
very closely indeed, but the signature may perhaps be 
sought in the statistics of the flow fields. 

I thank A. Louis, R. D. Groot and M. E. Gates for 
discussions. 
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FIG. 1. Interface density profiles for the three many-body 
DPD parameter sets in Table ^ The origin of the z-axis is 
arbitrary. In all cases there is a fairly sharp interface between 
a liquid of moderate density on the left, and an extremely 
dilute vapour phase on the right. 




FIG. 2. Density maps in bubble coarsening regime for set 
'6' in Table |, at reduced times t/Tn « 30 (left) and 120 
(right). Grey scale runs linearly between p = (black) and 
p > 6 (white). 
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FIG. 3. Scaling collapse of mean bubble area (A) as a func- 
tion of time, non-dimensionised by the hydrodynamic length 
and time scales, for the three many-body DPD parameter 
sets in Table ^ (upper data curves). Also shown is the corre- 
sponding reduced data for droplet area growth in binary liq- 
uid demixing, for the two conventional DPD parameter sets 
in Table ^ (lower data curves; the lowest being set '4b'). 
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Set 


A 


B 


rd 


P 


V 


0" 


Lh = rf Ipff 


Th = ??Vpo-^ 


'4' 
'5' 
'6' 


-40 
-40 
-40 


27 
18 
12 


0.75 
0.75 
0.75 


3.954(1) 
4.950(1) 
6.333(1) 


3.19(3) 
2.53(3) 
1.69(2) 


2.83(7) 
5.29(8) 
9.44(6) 


9.09 X 10"^ 
2.44 X 10"^ 
4.78 X 10"^ 


1.03 X 10" 
1.17 X 10-1 
8.55 X 10"^ 



TABLE I. Parameter sets for many-body DPD, for vapour-liquid coexistence, and properties of 

the liquid phase determined by simulation (in all cases, the figure in brackets is an estimate of the 
error in the final digit). The vapour phase density is vanishingly small. The final two columns show 
the hydrodynamic length and time scales. Units are DPD units, where m = Vc = kT = 1. 



Set All = A22 A12 p 77 a Lh = / po Th = "n'' j pa'^ 



'4a' 


5 


40 


4 


1.15(5) 


3.90(4) 


8.48 X 10"'' 


2.50 X 10-^ 


'4b' 


20 


50 


4 


1.10(3) 


6.55(4) 


4.62 X 10"^ 


7.76 X 10"^ 



TABLE II. Parameter sets for ordinary DPD, for symmetric binary liquid mixtures.. 
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